6 VACUUM IN A RANDOM MATRIX MODEL 
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Inspired by recent lattice calculations, we model certain aspects of the 6-vacuum 
using a matrix model with gaussian weights. The vacuum energy exhibits a cusp 
at d < TT that is sensitive to both the accuracy of the numerical analysis and the 
maximum density of winding modes present in a finite volume. 

Recent lattice study of the CP'^ model on the latticed have revealed an 
intriguing result: the energy as a function of the CP breaking 9 angle levels off 
at large 6. The result is very interesting, since in that model the string tension 
for small charges is related to the ^-derivative of the vacuum energy. Ajjlateau 
indicates the loss of confinement. The result was reexamined in u where 
the authors pointed out to some subtleties regarding the technical methods 
applied. 

In this talk we will address these issues using a random matrix model 
with the partition function 

„,'^.»,)>' « 

where the averaging is carried over the matrices W, W'' using the weight 
exp (— n/2TrVF^VF) and over the matrix size n and a, with weights 
exp (— x^/2x*V^) and exp (— cr^/2cr»F), respectively. Here is a complex 
asymmetric n+ x n_ matrix, n = n+ +n_, and a±x = '2n± — (n). The mean 
number of zero modes (n) is either fixed from the outside or evaluated using 
the gaussian measure. Here and for simplicity we use the quenched measure 
without the fermion determinant to fix (n). Throughout, we set the value of 
the quark condensate to S = 1 in the chiral limit. p. 

Eq. is borrowed from the effective instanton vacuum analysis U where 
n+ counts the number of right-handed zero modes , and n„ the number of left- 
handed zero modes. The number of exact topological zero modes is commen- 
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surate with the net winding number carried by the instantons and antiinstan- 
tons. By analogy with 0, x* ^-nd cr* will refer to the unquenched topological 
susceptibility and particle compressibihty, respectively, with a1 — 12n,/llAfc 
and ~ (n) /V the mean density of zero modes. If the compressibility ct» is 
assumed small in units of S = 1, then typically n ~ (n). 

quenched case: In the case of the quenched partition function {Nf — 0) 
one probes the nature of the gaussian measure. We start discussing the case 
where (n) = oo with no restriction on the value of n, and hence no restriction 
on the value of x- Then, we discuss the case where (n) is large but finite, so 
that Ixl < n with typically n ^ {n) for a peaked distribution in n. 

When the sum is unrestricted and infinite, using Poisson resummation 
formula we have 

+ 00 

ZQiO)= e-^^^*('-2"'=)' =03(^/2,6-^). (2) 
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with 6 being the third elliptic function and r — l/(2V^x*)- The result is 
manifestly 27r periodic. The vacuum energy, Fq{6) — ~IiiZq(6)/V as 1^ ^ oo 
is simply Fq{9) — minix* {6 + mod27r)^ in agreement with the saddle-point 
approximation. This result is in agreement with the result using large Nc 
arguments □, and recent duality arguments tl. We observe that the cusp at 
6 = TT (mod 27r) sets in for y = oo. 

In the matrix model being considered the sum over x is restricted to 
1x1 < N , with N — max n. We denote by n = N/V the maximum density of 
winding modes. While in general n 7^ n*, for a peaked distribution in n (small 
compressibility ti* ) we expect n n» . This will be assumed throughout unless 
indicated otherwise. Hence 

Approximating the sum in (||) by an integral and evaluating it by saddle 
point we obtain Zq ^ e^'^XtS /2 apparently in agreement with the infinite 
sum. However, the Euler-MacLaurin summation formula shows that deviation 
from this result in the case of a finite sum is expected for 9 > 0^ ^ n/x* u. 
Indeed, this is confirmed by detailed numerical calculations of the vacuum 
energy as shown in Fig ^ (left). For N — 250 and x* = the double precision 
(16 digit) numerics (circles) breaks away from the saddle point approximation 
(solid line) at 9/tt ~ 0.2 for both n = 1 and n = 4, while the high-precision 
(64 digits) calculations (dashed line) agree with the saddle point result at 
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n = 4 but break away at Qj-K ^ 0.3 « l/rTX* at n = 1. On the right we show 
the numerical result for the CP^ model lil showing a similar behavior. 
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Figure 1. Free energy in a quenched model (left) and in the CP^ model (right). See text. 

unquenched case: Analytical ancLpumerical studies of the unquenched 
matrix model for N f = 1 and N j ~ 2n lead to results similar to the ones 
observed in the quenched case. The mean- field, result breaks down for Q > 
Qc ~ n/xtop. For small quark masses Xto-p ~ m 13, the breakdown is observed 
only for ensembles with very small maximal winding number. 

In this talk we have shown that the vacuum energy of a chiral matrix 
model developes a cusp at < tt that is sensitive to the numerical accuracy. 
The cusp persits at high accuracy if the maximum winding number considered 
is low, in disagreement with mean- field results. Similar observations in current 
lattice simulations should therefore be interpreted with care. 
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